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Geometric-Phase-Effect Tunnel-Splitting Oscillations in Single- 
Molecule Magnets with Fourth-Order Anisotropy Induced by 
Orthorhombic Distortion 



M. S. Foss-Feig'''' and Jonathan R. Friedman*''' 
Department of Physics, Amherst College, Amherst, MA 01002-5000 



PACS 7 5 . 4 5 . + j - Macroscopic quantum phenomena in magnetic systems 
PACS 7 5 . 5 . Xx - Molecular magnets 

PACS 3 . 6 5 . Vf - Phases: geometric; dynamic or topological 

Abstract. - We analyze the interference between tunneling paths that occurs for a spin system with both 
fourth-order and second-order transverse anisotropy. Using an instanton approach, we find that as the strength 
of the second-order transverse anisotropy is increased, the tunnel splitting is modulated, with zeros occurring 
periodically. This effect results from the interference of four tunneling paths connecting easy-axis spin orien- 
tations and occurs in the absence of any magnetic field. 



Geometric-phase effects play an important role in spin dy- 
namics. Notably, a geometric phase lies at the heart of spin- 
parity effects, such as Kramers degeneracy [1-3]. Tunneling 
of a spin (or magnetic particle) between degenerate orienta- 
tions can be modulated by such geometric-phase effects via the 
interference between multiple tunneling paths, with complete 
suppression (or "quenching") of the tunnel splitting occurring 
when tunneling paths interfere destructively. 

Garg has considered the case in which the tunneling of a spin 
with biaxial anisotropy can be modulated by the application 
of a magnetic field along the hard axis [4,5]. Specifically, he 
studied the Hamiltonian 



D 



(1) 



for spin J and anisotropy constants D > E > and = 
gfisHx- He showed periodic quenching of the tunnel split- 
ting as a function of magnetic field due to geometric-phase 
interference between tunneling paths. Several extensions and 
variants of this problem have been studied theoretically [5-18]. 
In a ground-breaking experiment, Wernsdorfer and Sessoli [19] 
found clear evidence of this effect in the single-molecule mag- 
net Feg, which is reasonably well described by Eq. 1. 

The case originally studied by Garg involved no magnetic 
field along the longitudinal (z) direction. Tunneling can occur 
between the metastable ground state | J) and an excited state 
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I — J + n) when the applied field brings them into resonance, 
i.e. when = giJ-sHz — nD, n = 0, 1, 2, . . . , 2 J — 1 [20]. 
In Feg, it was discovered that at these resonance longitudi- 
nal fields, a transverse field along the hard axis continued to 
produce interference effects [19]. Thus, at certain values of 
hz and h^, destructive interference occurs, suppressing tunnel- 
ing. These so-called diabolical points in hz-h^ parameter space 
have been studied theoretically by several authors using differ- 
ent approaches [5-7, 10, 1 1, 15, 17]. 

The Mni2 single-molecule magnet has a four-fold transverse 
magnetic anisotropy and, like Feg, displays resonant tunnel- 
ing between two easy-axis orientations [20]. One variant of 
this molecule, Mni2-tBuAc, appears to display the four-fold 
symmetry with high accuracy [21-23]. Here we consider the 
interference that occurs between tunneling paths in such a sys- 
tem and how that interference can be modified by the presence 
of a second-order transverse anisotropy perturbation. Such a 
perturbation could potentially be induced by the application of 
uniaxial pressure to a sample of Mni2. We find that the tunnel 
splitting is periodically quenched as a function of the strength 
of the perturbation. This interference effect takes place in the 
absence of any magnetic field. In addition, we find a pattern of 
diabolical points in parameter space when a longitudinal field 
hz is applied. 

Several investigations have studied the inclusion of a fourth- 
order transverse anisotropy as a perturbation in Eq. 1 that af- 
fects the number and positions of quenches produced by the 
external field [5, 11, 16, 17, 19]. Other work has studied how 
a transverse magnetic field modulates the interference among 
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the four tunnehng paths defined by a fourth-order transverse 
anisotropy [12, 13]. In the present study the fourth-order 
anisotropy term is the primary transverse anisotropy in the 
problem, producing four interfering paths, but the interference 
is modulated by the strength of the second-order anisotropy, 
leading to periodic quenches as that term is varied. In particu- 
lar, we consider a spin governed by the Hamiltonian 

H = d(^J^ - jfj +k(.jX + Jl) + A ( J2 - jfj , (2) 

where the J's are standard spin operators. We restrict the val- 
ues of k and A such that the z axis is the spin's easy axis (see 
Appendix). 

We calculate the ground-state tunnel splitting using an in- 
stanton approach, valid in the semi-classical limit of large 
J. Since we are considering tunneling between energy min- 
ima, there can be no classical trajectories connecting the end- 
points of the motion. The quantum trajectories with the 
largest amplitude are the instantons, obtained from the Euclid- 
ian (imaginary-time) Lagrangian 

Le = L{t^-iT) = ij {cose -1) ~n{e,(t>), (3) 

dr 

where 7i (61, 0), the Hamiltonian expressed in terms of spherical 
angles, is identified with the spin's classical energy Eci{0,(j)) 
defined in Eq. 6 below. Since the instanton paths are extrema 
of the action in imaginary time, the total energy must be a con- 
served quantity along each path. 

The ground-state (zero-temperature) tunneling rate is given 

by 

A = ^7,e-^--, (4) 

i 

where the sum is over all instantons and the pre-factors ji arise 
from integrating the fluctuations around the instanton paths 
[24]. The Euclidean action Se for an instanton path is 

Se = - LecLt 

= ij {i-cose)d(j)+ H{e,4>)dT (5) 

The first term has the form of a Berry phase: The difference 
in the imaginary part of the action between two paths is pro- 
portional to the real part of the solid angle (on the complexified 
unit sphere) circumscribed by the two paths. Since the last term 
is the total energy, it is a constant for an instanton path and, by 
construction, zero for the Hamiltonian, Eq. (2). 

For the case of Eq. (2), the spin dynamics can be under- 
stood in terms of four instanton paths along the complexified 
unit sphere. When A = 0, the spin's energy has four saddle 
points located in the x-y plane, with instanton paths connect- 
ing the +z and -z directions virtually "passing through" these 
saddle points, as shown in Fig. la, where the solid red curves 
show numerically calculated instanton paths with only the real 
parts of 9 and </> displayed. For A > 0, the energy landscape is 
changed, as shown in Fig. lb, and the instanton paths are pulled 



towards the y axis. Reflection symmetry in the x-z and y-z 
planes requires that all four paths have the same amplitude, but 
each has a different geometric phase that leads to interference. 
We therefore concern ourselves with calculating the imaginary 
part of the action. The interference effect can be understood by 
considering, say, the top two paths in Fig. 1 . The solid angle be- 
tween these paths decreases with increasing A, modulating the 
interference between their tunneling amplitudes. An identical 
effect occurs for the bottom two paths. 

For ease of computation we measure energy in units of D 
(i.e. we set D = 1). We choose slightly unconventional spher- 
ical coordinates, with the polar angle 9 measured from the x 
axis and the azimuthal angle </> measured from the y axis to 
avoid having coordinate singularities at the poles. In terms of 
these angles, the classical energy associated with Eq. (2) is 

E^i{9,(j)) = J2(l - sin^eisin^^) 

+\j'^{cos^9-iim^9cos'^ (f) 
+2kJ^{cos^ 6 + sin^ 9 coa'^ </> 
-Bcos^ 61 sin^ 61 008^0), (6) 

which has two degenerate absolute minima (with Eci = 0) at 
6' = 7r/2,0 = ±7r/2. 




Fig. 1: (Color online) Three-dimensional plot of energy landscape for 
a spin governed by Eq. 6 with a = kj^ =0.1 and A = (a) or 
A = 0.2 (b). The energy minima are along the z axis (into and out 
of the page) while the x and y axes are the hard axes. The solid red 
curves show the four instanton paths with only the real parts of 6 and (j> 
displayed. These paths were numerically calculated using Eqs. 7 and 
10. A rather large value of a is used here to ease visual representation. 
ForMni2, a 3 x 10"^. 

We can quickly gain a large amount of information about 
the instantons by solving Eq. (6) for the 9{(f)) paths that con- 
serve energy. Making the substitution p — J cos 9, the conju- 
gate momentum to (p, the instantons must satisfy the constraint 
Eel — 0, which has four solutions: 

P±,±{'P)=±J^ , (7) 

where the subscripts refer to the choices of signs and 
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A{(f>) = 4a cos^ 4> + 12a cos^ (j) - sin^ (j) ~ X cos^ - A 
5(0) = -(l28a2 + 16a(A-2) + (A-l)2)cos'*0 
+2 (1 - A^ + 8a(A + 2)) cos^ - (A + 1)^ 
(7(0) = 40(003"* + 6cos2 0+1), (8) 

with a = kJ"^. 

Two of the solutions in Eq. (7), p±.+, do not have values of 
p = at the endpoints = ±7r/2. While such "boundary- 
jump instantons" [11] are allowed solutions, numerical evalu- 
ation of these solutions shows that their action is much larger 
than for the other two solutions and therefore they give negli- 
gible contribution to the tunnel splitting. We neglect them in 
further discussion. 

The other solutions p±,- for < a < 1/2 and |A| < Ac = 
Ay/a — 2a^ have four branch points along the real-0 axis at 
— ±7r/2 and at the zeros of -6(0). Evaluating the action 
along the real-0 axis is then potentially problematic and the in- 
tegration contour must be chosen with some care. A judicious 
choice of branch cuts is shown in Fig. 2. We find the instanton 
trajectories by considering the (imaginary) time dependence of 
the coordinates. 



Im (j)/7i 
0.15 r 




-0.15 L 



Fig. 2: Instanton contour in the complex-0 plane for a = 10~* and 
A = 0. Branch points and branch cuts are shown. The contour changes 
only slightly as a is varied by orders of magnitude. 

The equations of motion for the variables and p can be 
found using Hamilton's equations: 

= -2i,/p(4a/ + J^sin^0 + AJ^(1 + 008^0) + 

12a(2p2 - J2) cos^ + 4a(j2 - p^) cos^ 4))/j'^ (9) 

p = -isin(20)(j2-p2)(2a(j2-7p2) + J^(l-A) + 
2a(j2 -p2^cos(20))/j2. (10) 

Eq. (9) can be decoupled by substituting one of the solutions 
p±,- for p. The equation can then be solved numerically to 
find the real and imaginary parts of 0(t). For each choice of 
p±- we find two solutions, one with endpoints at = ±7r/2 
(passing through = 0) and the other with endpoints = 7r/2 



and 37r/2 (passing through = tt) Four more solutions can 
be found by noting that Eqs. (9) and (10) are invariant under the 
transformation (0 ^ — 0*,p ^ p*) followed by complex con- 
jugation of the equations. This yields a total of eight solutions, 
four of which are instantons and four anti-instantons. These 
can be sorted by using Eq. 9 to determine the sign of at the 
midpoint of each solution. Table 1 indicates the Euclidean ac- 
tion integrals for the four instantons, with the limits designating 
the direction of the path about the polar axis. 



Instanton 


Euclidean Action 


1 




2 




3 




4 





Table 1: Instanton Actions: the upper bound — 7r/2 versus 37r/2 is 
intended to convey whether the instanton runs through or vr. 

By symmetry, only one of these integrals need actually be 
evaluated. The contour in the complex-0 plane associated with 
one of these solutions (Instanton 1) is shown in Fig. 2. That 
the contour begins and ends at a branch point poses no compli- 
cations since the integrand vanishes smoothly at those points. 
With the branch cuts shown, the contour can be continuously 
deformed onto (or within e of) the real-0 axis. Along this axis, 
the real part of p± _ (the part we are interested in) is only non- 
zero between the inner two branch points. The imaginary part 
of the Euclidean action (for Instanton 1) can then be written as 

Si = -J{tt + I) 
^ = jI Re[p+,_(0)]d0, (11) 

where ±6 are the positions of the inner branch points. Proceed- 
ing similarly for the other instantons, the total tunnel splitting 
can be calculated using Eq. (4): 

A(A) oc e-^«[e-"^"(e-''^^ + e*'^^) 

= 46^^^^ cos ( Jtt) cos ( JI) , (12) 

where Su is the real part of the Euclidean action and the same 
for all instantons. The factor cos (Jtt) ensures that the tun- 
nel splitting is zero for half-integer spin, in accordance with 
Kramers' theorem [1,2]. It arises from the interference be- 
tween diametrically opposite paths (Fig. 1) that remains con- 
structive (destructive) for integer (half-integer) values of J as 
A is varied. The cos (JI) factor is of chief concern here. Since 
the integral / depends on anisotropy parameters a and A, varia- 
tions in these parameters can modulate the tunnel splitting, with 
complete suppression of tunneling when JI = {2n + 1)tt/2 

' While our numerical calculations show that each is a solution, we have not 
proven that these are the only solutions. However, we do not expect more than 
four instanton solutions with equal magnitude actions given the symmetry of 
the problem. 
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for integer n. For fixed a, as A is increased from zero, the in- 
ner branch points move towards each other and the value of / 
decreases, resulting in repeated suppression of tunneling until 
A = Ac = 4\/a — 2a^ at which point 6 = and the value of / 
goes to zero. 

An approximate expression for / can be obtained by treating 
the integrand as a function of two independent parameters (a 
and A/ Ac) and expanding in a power series in a. The resulting 
integrals can be done by elementary means. We obtain 



/=^(l-^).(a,A), 



(13) 



where g (a, A), up to order a^, is 
5 (a, A) 



1+ 11+ 
Ac \ Ac 




(14) 



Fig. 3 compares the results of Eqs. (13) and (14) (solid blue 
curve) to exact numerical evaluation of / (points) for a — 



2-V2 



0.073, the largest value of a for which I can be mean- 



ingfully evaluated for all A < Ac (see Appendix). Even at this 
relatively large value for a, Eq. (13) is clearly an excellent ap- 
proximate expression for /, although for such large values of a 
and A the instanton approximation is questionable since the bar- 
rier has nearly disappeared. The A-dependence of / is clearly 

dominated by the ^1 — factor in Eq. (13), responsible for 
the nearly linear behavior seen in the figure. To illustrate the 
geometric interference effect, the figure also shows the behav- 
ior of the magnitude of the cos (J/) factor in A as a function of 
A (red dashed curve) for J — 10. We note that when J is an odd 
integer, the zeros in A start at A = (because / = 7r/2 at this 
point), as expected by general symmetry considerations [2]. 
Since the argument of cos {J I) ranges from zero to Jtt/2, the 
number of zeros is [( J + 1)/2J . 

We next investigated the effect of a longitudinal field on the 
geometric -phase interference. Such a field introduces a term 
—hzSz into the Hamiltonian, breaking the symmetry between 
up and down directions of spin (i.e. raising the energy of one 
minimum with respect to the other) but preserving other ro- 
tational and reflection symmetries. In particular, the four in- 
stanton paths maintain the same amplitude and will continue to 
interfere. 

Similar to what was found for the biaxial spin, we deter- 
mined that a spin described by Eq. 2 has an array of diabol- 
ical points. At certain values of longitudinal field, varying A 
causes the tunnel splitting to oscillate, with perfect degenera- 
cies occurring at a discrete set of points in X — space. The 
effect can be most readily addressed by invoking the argument 
of Bruno [17] that tunneling between \J) and \ — J + n) can 
be mapped approximately onto the problem of tunneling be- 
tween the ground states | J) and | — J) of an effective spin 
J = J — n/2 in zero magnetic field. Making the substitution 
J ^ J in Eqs. (12) and (13) and keeping in mind that both a 
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Fig. 3: (Color online) I/tt (blue, solid curve and points) and | cos (J/) | 
(red, dashed) as a function of A/ Ac for a = '^^^ and J — 10. The 
points are the result of numerical evaluation of / whereas the solid blue 
curve is obtained using Eq. (13). The value of a used is the largest one 
for which the problem is well defined for all values of A < Ac (see 
Appendix). 

and Ac are functions of J, we can generate the full set of diabol- 
ical points. If a is small, we can approximate g {a, A) sa 1 and 
get the following simple expressions for the diabolical points. 



nD, 



A = AdkD-2{kJY 



J -21-1 



(15) 



where I = 0, 1, ... , [(J — 1)/2J and n is an even (odd) inte- 
ger when J is a (half-)integer. For other integer values of n, 
Eq. (12) yields exactly zero tunnel splitting for all values of A, 
as expected in the absence of a transverse field. In these ex- 
pressions, we have restored the anisotropy parameter D. The 
first expression is equivalent to the condition for resonant tun- 
neling [20] in which the ground state associated with the up di- 
rection is resonant with the nth excited state in the down direc- 
tion. Fig. 4 shows a comparison between the results of Eq. (15) 
and exact numerical diagonalization of the Hamiltonian using 
J = 10 and k/D = 2 x 10"^, a value close to that of Mni2. 
The agreement is obviously very good. 

We note that the pattern of diabolical points in Fig. 4 resem- 
bles the pattern found in the biaxial system [5-7, 10, 1 1, 15, 17], 
where the diabolical points lie on a perfect lattice in hj. — 
space. In contrast, the diabolical points in Fig. 4 lie in the X-h^ 
plane and only form a perfect lattice in the limit of fc — > 0. 

Finally, we consider the feasibility of observing the pro- 
posed effect. For Mni2Ac, a « 3 x 10"'^, which gives 
Ac w 0.2D — 0.11 K. Density functional calculations that in- 
clude molecules of solvation in a trans configuration estimate 
that a 3% strain along the axis defined by the solvent molecules 
can produce changes in A as large as 19 mK [25]. To the best 
of our knowledge, the elastic modulus of Mni2Ac has not been 
measured, but it is unlikely that it could be much larger than 
100 kbar One might then need uniaxial pressure as high as 17 
kbar to reach Ac, although the interference effect would become 
manifest at substantially lower values of A. Several experimen- 
tal studies have investigated the effect of hydrostatic pressure 
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move towards each other (see Fig. 1). They merge together 
when A = 8a (dashed line in Fig. 5). Remarkably, even when 
there are only two saddle points in the potential landscape, there 
remain four distinct instanton paths. The solid angle between 
pairs of paths vanishes when X — Xc (solid curve), which does 
not correspond to any notable change in the energy landscape. 
The vertical dashed line represents the values in this parameter 
space used for the results presented in Fig. 3, i.e. a is such that 
Xc = 1. 

Finally, when A > 1 each energy minimum at a pole bifur- 
cates into two minima, one tilted towards the +y direction and 
the other towards the -y direction. This substantially changes 
the nature of the problem. The analysis described in this paper 
applies only to the unshaded region in Fig. 5 and would need 
modification to handle the energy landscapes corresponding to 
the shaded regions. 



Fig. 4: Diabolical points in X ~ hz parameter space with k = 
2xlO~^DandJ = 10 (i.e. Mni2) calculated using Eq. (15) (squares) 
and by exact diagonalization of the Hamiltonian ( x 's). Points at neg- 
ative values of hz and A were obtained by reflection in the horizontal 
and/or vertical axis. 

on Mni2 [26-29] and have found significant changes in D with 
pressures as high as 12 kbar. None of these investigations in- 
volved uniaxial pressure, which would be needed to observe the 
interference effect predicted herein. We note that the ground- 
state tunnel splitting for Mni2Ac is extremely small. However, 
the splitting increases exponentially for progressively higher 
levels, suggesting that a pressure-induced effect would be ob- 
servable in thermally excited states, i.e. by studying the mag- 
netic relaxation in the thermally assisted tunneling regime [20] 
as a function of pressure. Discussion of the geometric -phase 
effect in excited levels will be presented elsewhere. 

In summary, we predict a geometric-phase interference ef- 
fect involving four instanton paths that can be modulated by 
the strength of the second-order transverse anisotropy, leading 
to a nearly periodic suppression of tunneling in the absence of 
any external field. In principle, this effect can be experimen- 
tally realized through the application of uniaxial stress along 
one of the hard-axis directions of a four-fold symmetric single- 
molecule magnet, such as Mni2. 

Appendix. - Here we discuss how the energy landscape 
for a spin governed by Eq. 6 depends on the anisotropy param- 
eters a and A. (As in the main text, we measure these anisotropy 
parameters in units of D and define a = kJ^.) Using standard 
techniques from multivariable calculus we determined the min- 
ima, maxima and saddle points of the landscape. 

For small values of a and A (unshaded region in Fig. 5), the 
landscape resembles that shown in Fig. 1 with global minima 
at the poles (±z directions) and saddle points along the equator 
(x — y plane). As a increases, the saddle points become deeper 
and for a > 1/4 the saddle points transform into local minima 
on the equator (right shaded region in Fig. 5). This limit is far 
from relevant for Mni2, where a « 3 x 10^"^. 

Increasing A causes pairs of saddle points on the equator to 
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Fig. 5: Parameter-space plot indicating the essential features of the 
spin's energy landscape. The unshaded region represents the locus of 
parameters to which the analysis in this paper applies. In the shaded 
region on the right, the saddle points on the equator have become local 
minima. In the shaded region above A = 1, the minima at the poles 
have bifurcated. The dashed line indicates where equatorial saddle 
points merge pairwise. The solid curve shows the behavior of Ac and 
is included for reference. 
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